www.gusucode.com > 基于matlab编程蒙特卡罗仿真设计源码程序 > 基于matlab编程蒙特卡罗仿真设计源码程序/第四讲 随机数的生成及随机变量抽样/liti9.m

    mm=10000
Randnum1=unifrnd(-1,1,1,2*mm);
Randnum2=unifrnd(-1,1,1,2*mm);
xRandnum=zeros(1,mm);
yRandnum=zeros(1,mm);
s=Randnum1.^2+Randnum2.^2
ii=1;jj=1;
while ii<mm
    if  s(1,jj)<=1;
       xRandnum(1,ii)=Randnum1(1,jj);
       yRandnum(1,ii)=Randnum2(1,jj);
       ii=ii+1;
    end
    jj=jj+1;
end
[Y,JJ] =sort(xRandnum);
plot(xRandnum(JJ),yRandnum(JJ),'.')